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Abstract: We present a supersymmetric model in which the observed fermion 
masses and mixings are generated by localizing the three generations of matter and 
the two Higgs fields at different locations in a compact extra dimension. Super- 
symmetry is broken by the shining method and the breaking is communicated to 
standard model fields via gaugino mediation. Quark masses, CKM mixing angles 
and the \x term are generated with all dimensionless couplings of 0(1). All dimen- 
sionful parameters are of order the five- dimensional Planck scale except for the size 
of the extra dimension which is of order the GUT scale. The superpartner spectrum 
is highly predictive and is found to have a neutralino LSP over a wide range of pa- 
rameter space. The resulting phenomenology and interesting extensions of the model 
are briefly discussed. 
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1. Introduction 



Particle physics is littered with energy scales. The known particles have masses 
which are spread over numerous orders of magnitude. Most hadron masses hover 
near the scale of presumed quark confinement. However, the masses of the eight 
light mesons or pseudo-Goldstone bosons, may be parameterized by explicit chiral 
symmetry breaking in the fundamental theory. In QCD, this breaking for the most 
part is due to small quark masses. 

In the standard model (SM), non-zero quark and charged lepton masses require 
the electroweak symmetry to be broken. This breaking (EWSB) is accomplished by 
allowing a scalar field (Higgs boson) to have a vacuum expectation value (VEV) thus 
giving masses to the W and Z bosons, the gauge fields of the electroweak interaction 
(minus the photon). One might expect the masses of the charged fermions to be of 
similar magnitude as the scale of electroweak symmetry breaking. However, instead 
the masses extend more than five orders of magnitude below the weak scale. One 
(substantially explored) explanation for this large hierarchy of masses is the existence 
of some symmetry broken at high energies producing masses suppressed relative to 
the weak scale. 

A more serious hierarchy problem in the SM is the instability of the Higgs mass, 
and thus the electroweak scale, with respect to quantum corrections. A theory with 
supersymmetry (SUSY) softly broken at the weak scale, such as the minimal super- 
symmetric standard model (MSSM), can stabilize the Higgs mass because in such 
theories all quadratically divergent quantum corrections vanish. However, generically 
these theories contain contributions to flavor changing neutral currents (FCNC) and 
CP violation which far exceed current experimental bounds 0. Again, symmetries 
could restrict flavor violations among supersymmetry breaking terms. For instance 
flavor symmetries which also restrict the soft terms , gauge symmetries which me- 
diate supersymmetry breaking |J, or some combination which results in partially 
aligned heavy first two generations |4], [|. 

Recently, Randall and Sundrum have suggested a new way to explain small 
couplings without appealing to symmetries. They were able to forbid generally fla- 
vor violating non-renormalizable operators that mix MSSM fields with fields in a 
supersymmetry breaking sector by spatially separating the two sectors in a small 
extra dimension. All soft terms appear due to contributions coming from the su- 
perconformal anomaly || [7|. While the soft terms are sufficiently flavor diagonal in 
the minimal scenario, the sleptons are tachyonic, and thus break the electromagnetic 
symmetry. Various model-building scenarios have appeared in the literature which 
attempt to fix this problem ||. 

Arkani-Hamed and Schmaltz have since shown that localizing fields in extra 
dimensions at distances of order unity with respect to the fundamental scale can 
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easily produce exponentially small Yukawa couplings and at the same time suppress 
proton decay due to the small overlap the their wave functions ||. This mechanism 
is especially useful in models in which large compact extra dimensions solve the 
hierarchy problem by bringing the fundamental Planck scale down to just above the 
weak scale [pTOp . In general, this method presents an interesting alternative to the 
usual spontaneous flavor symmetry breaking scenarios JI2 . 



In this note, we present a model for fermion masses where Yukawa couplings 
are suppressed due to the localization of fields in an extra spatial dimension. Our 
model differs significantly from the Arkani-Hamed/Schmaltz model in three ways*: 
(i) our extra dimension is small (tens of Planck lengths) and supersymmetry solves 
the hierarchy problem, thus avoiding bounds from flavor- violating neutral current in- 
teractions induced by the relatively light Kaluza-Klein (KK) excitations of the gauge 



bosons | (ii) Yukawa couplings are small due to the position of each generation 
relative to the (localized) Higgs fields and not due to the different splittings of left 
and right handed fermions [[OJ, and (iii) our gauge fields fill the entire space in the 



extra dimension thus relieving us of the difficult task of localizing gauge fields in field 
theory [14], [15]. 



Because our gauge fields live in the bulk, localized supersymmetry breaking 



produces the conditions necessary for recently proposed []I8[ [19] gaugino-mediated 
supersymmetry breaking (gMSB), arguably the simplest way to mediate supersym- 
metry breaking while avoiding all phenomenological flavor constraints. In addition, 
by using a particular variation of the "shining" mechanism of Arkani-Hamed, et. al. 
[fZOfl , we localize supersymmetry breaking close to the Higgs fields making it possible 
to both generate a /x term and insure that the right-handed stau is not the lightest 
super symmetric particle (LSP). 

Our fields are localized by generation and thus are consistent with a supersym- 
metric SU(5) grand unified theory (GUT). This has implications for charged lepton 
masses and thus we shall assume that GUT-breaking is responsible for producing the 



correct leptonic spectrum [[H]] . A small extra dimension opens new possibilities with 
regards to GUT models, but we leave this and the leptons for future work. 

The paper is organized as follows. Section 2 describes all of the elements nec- 
essary for a successful flavor model in this context. We find the model to be quite 
constrained and predictive. Section 3 describes the incorporation of supersymmetry 
and supersymmetry breaking, including a brief review of gMSB and some attractive 
modifications with respect to the fi term. Section 4 describes the phenomenology of 



*Dvali and Shifman jl6| also considered localizing complete generations and Higgs fields in an 
extra dimension and suggested a way supersymmetry breaking could be included via non-BPS 
brane configurations. Also, Gherghetta and Pomarol have recently suggested that small Yukawa 
couplings could arise from localizing fermions at different positions in a slice of anti-de Sitter space 
and discuss variations of supersymmetry breaking in such scenarios ]i~7fl . 
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the model with respect to the quarks and CKM matrix as well as the supersymmet- 
ric spectrum. Section 5 discusses possible future enhancements of this model and an 
appendix describes the localization of an N = 1 chiral multiplet . 



2. Flavor from Extra Dimensional Overlaps 

In this section we show how, with a small extra dimension, to generate small 
fermion masses from 0(1) Yukawa couplings. The fermions are localized with respect 
to the extra dimension and their separations from each other and the Higgs fields are 
also 0(1) (in Planck units). Supersymmetry does not play a role in this discussion, 
other than to motivate the existence of two Higgs doublets, and thus the results of 
this section may be applied generically to any supersymmetric theory. The details 
concerning the localization of a zero mode of a chiral superfield with a Gaussian 
profile are presented in Appendix In this work we will only consider the quark 
masses and mixings. The lepton masses may be obtained from a straight-forward 
generalization of these results. 

2.1 Gaussian Localized MSSM fields 

We take as our starting point a five- dimensional (5d) GUT. The fifth dimension is 
small, with associated mass scale M c of order 1/100 times the four-dimensional (4d) 
Planck scale. Thus, at energies below M c , there is a 4d effective theory description 
of the resulting dynamics. The fermion masses arise from 5d superpotential terms, 
written conveniently in the language of 4d iV = 1 supersymmetry [[2(J as, 

J dy J d 2 6j2 I ^H u Q t U^ + ^H^dA + H.c., (2.1) 



where H u and are the chiral superfields containing the up- and down-type Hig- 
gses, Q is the quark SU(2) doublet, and U c and D c are the up- and down-type 
quark SU(2) singlets. The 5d Planck scale M* is related to the 4d Planck scale 
by M* = (Mp/L) 1 / 3 , y is the coordinate parameterizing the compact dimension and 
powers of M* have been inserted such that are dimensionless. This superpotential 
violates the N — 1 supersymmetry in 5 dimensions but is invariant under half of the 
supersymmetry transformations which correspond to N = 1 supersymmetry in 4 di- 
mensions. One could imagine explicit breaking in a microscopic theory where at least 
some of the fields (either quarks or Higgses) in ( |2.1|) live on "3-branes" which break 
translation symmetry in the 5th coordinate along with half of the supersymmetries^. 

^This assumption could lead to fields with a simple exponential fall-off rather than the Gaussian 
profile described below. We have checked to see that successful models could be produced with 
simple exponentials and found similar results to those described here. Thus, for simplicity, we 
assume all fields have Gaussian wave functions. 
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We also assume it is possible to localize only a single left-handed zero-mode, without 
also introducing a right-handed one. We will not explore this subtle but difficult issue 
and assume it can be accomplished. Localizing a field using the method outlined in 
the appendix also explicitly breaks half the supersymmetries and again we take this 
as a requirement of the form of explicit "N — 2" supersymmetry breaking and leave 
this issue for future work*. In addition, we assume the chiral superfield component of 
the 5d gauge multiplet can be given an explicit mass on a 3-brane, or removed from 
the low energy theory in some way (for example, see [2q|). Here we take a bottom- up 
approach in constructing the model, allowing successful phenomenology to motivate 
the high-energy theory. 

The Higgs and fermion fields are localized in the fifth dimension with Gaussian 
profiles. We will assume that the zero modes of the quarks of a given family are fixed 
to same location in y, and that the Gaussian profiles of the zero modes for all three 
families have a common width 2/( 2 (which we take to be about the Planck scale). 
This is consistent with GUT scenarios and would arise if all matter was localized 
by the same mass function (or VEV) whose slope doesn't significantly vary over the 
positions of the localized fields. The localization of each family around a different 
point in y can be accomplished by giving each family's hypermultiplets different 
constant mass terms in addition to the single mass profile which results in the zero 
modes || . Thus, provided the extent of the extra dimension is large enough that the 
deviations from this Gaussian profile are small, the zero mode profile for fermion j 
is given by, 

m= f^Ve-^) 2 . (2.2) 



3 V n / 

The Higgs superfields are also localized in y, and we further allow them to have 
different widths from the fermions and from each other. For now we will assume 
that all of the Yukawa interactions Yij are exactly 1. We will see below that a phase 
(required by CP violation) will also be important in obtaining the correct mixing 
angles. 

The resulting low energy effective theory has exponentially suppressed Yukawa 
interactions that result from the overlap of any two fermion wave function with a 
Higgs wave function. We define our coordinate system such that the up-type Higgs is 
at y — 0. We measure distances in y in units of and thus the model is completely 
specified by the locations of the three families and down- type Higgs, h, I2, h, and lh, 
and the relative widths of the Higgses, r u = Ch u /C and = (,H d K- In terms of these 
quantities, the resulting 4d Yukawa interactions for up-type quarks with H u is, 



2 + r 



' {{l+rl){ll + l])-2l l l j ) 



;2.3) 



*We thank Martin Schmaltz for discussions on these issues. 
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where Y~ has been taken to be one, and £ to be M*. There are also interactions 
between the down-type quarks and H d of the same form, but with r u — > r d and 
h — > (h — h)- When the up- and down- type Higgs scalars obtain VEV's v u and v d , 
these interactions will provide Dirac mass matrices for the quarks. These matrices 
may be diagonalized by separately rotating the right- and left-chiral fields, resulting 
in six real quark masses, and the three mixing angles and one complex phase that 
make up the Cabibbo-Kobayashi-Maskawa (CKM) matrix. Thus, our first task is 
to see if there exists a suitable arrangement of parameters to match the low energy 
data. 

2.2 Fermion Masses and the CKM Matrix 

We will now outline a method to determine our flavor parameters, Zi, I2, h, lh, 
r u , and r d , in order to fit the low energy data. We start with the MS quark masses 



and (90% CL.) three generation CKM matrix |22 



m u (2 GeV) = 1.5-5 MeV, m c (m c ) = 1.1 - 1.4 GeV, (2.4) 
m d (2 GeV) =3-9 MeV, m b (m b ) = 4.1 - 4.4 GeV, 

m s (2 GeV) = 60 - 170 MeV, m t (m t ) = 161 - 171 GeV, 



/ 0.9742 - 0.9757 0.219 - 0.226 0.002 - 0.005 \ 
\V C km\ = 0.219 - 0.225 0.9734- 0.9749 0.037- 0.043 

y 0.004- 0.014 0.035 - 0.043 0.9990 - 0.9993 ) 

We will work at the top mass scale§ (using the three-loop QCD and one-loop QED 



renormalization group scaling factors of [23[] to find the light quark masses at m t 



The measured Z boson mass requires v = a/ v\ + v d = 246 GeV, while the ratio 
tan/3 = v u /v d remains unfixed by EWSB. We will treat tan/5 as a prediction of 
our model, and characterize the allowed range of tan f3 by what results from the set 
of model parameters which accurately predict the quark masses and mixings. The 
correct top mass in eq. (|2]3|) is obtained by fixing the magnitude of l 3 once r u is chosen. 
(One can think that ^(r^) is fixed by m t to be a function of a chosen r u ). In the 
spirit of our work, we want to invoke numbers of 0(1), and thus we allow the widths 
to vary at most between 1/2 and 2. For v u ~ v and r u ~ 1.5, this requires / 3 ~ 0.3, 
though there is some freedom to vary Z 3 and r u together. Requiring that the ratio 
m t /m c come out correctly requires us to choose an appropriate l2(h, r u) = hij'u)- 
For the particular numbers discussed above the working choice is I2 ~ 2.3. 



§The 5d Yukawa interactions in Section 2.1 are actually those at M c . However, as the quark 
masses evolve together between m t and M c , this can simply be corrected by rescaling all of the 5d 
Yukawa interactions. 
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Turning to the down-type quarks, the ratio m&/m s , combined with the already 
"chosen" value of r u and the "determined" values of Z3 and I2 fixes a combination of 
and lh- The magnitude of m& determines tan/3 in terms of the above parameters. We 
may now choose l\ in order to arrive at the correct CKM element V us . As it stands, 
we have chosen two parameters (r u and r^), and used 4 pieces of experimental data 
(m t , m t /m c , mb/m s , and V us ) to determine the remaining parameters li, Z 2 , Z 3 , and lh- 
It remains to determine whether we can accommodate the remaining experimental 
data: V u b, V c b, and the first family quark masses by varying and r u independently 
over the "reasonable" range of 1/2 to 2. (One can think that two of the remaining 
observables fix and r u , and the last two are predictions of the model). 

As it turns out, the answer is no. While we easily realize the correct order of 
magnitude for the remaining predictions, we cannot quite reach the experimental 
values. V c b is always at least a factor of two or so smaller than its measured value. 
This is not really a serious failing; we have set out to construct a model in which 
all of the 5d Yukawa interactions were 0(1), but we have tried to realize the model 
by taking the couplings to be strictly one. This will never produce the CP violation 
observed in nature [22], and thus it is obviously too naive. In fact, taking a simple 
ansatz that the 5d Yukawa interactions have the phase structure, 



(2.5) 



we find that it is quite easy to fit the low energy data for a range of lh and r^. 
This particular choice produces a CP violating phase of 0(1), in accordance with 
measured CP violation in the Kaon system. We have verified that the addition of 
phases to the 5d Yukawa interactions greatly extends the workable range of li, Z 2 , l 3 , 
lh, r u and r^. Thus, in our context the phases are not only critical for CP violation, 
but also to get the right magnitude for the CKM elements. We stress that the ansatz 
(|2.5|) is only one of a wide range of workable solutions, which is a general indication 
of the robustness of the result. 

As a particular example of a working parameter set, consider the model defined 
by r u = 1.5, r d = 0.67, li = 3.05, Z 2 = 2.29, l 3 = 0.36, and lh = —1.7. The Gaussian 
profiles for the zero modes are shown graphically in fig. [l]. The physical origin of 
this solution to the flavor problem is evident from the positions and widths of the 
profiles. The third family sits close to the up-type Higgs to produce the large value 
of mt- The first and second generations are rather close together compared to the 
third generation, in order to correctly generate V us > V c b > V u b- The different widths 
and locations of the Higgs profiles account for the fact that m t > mb and m c > m s , 
whereas > m u . The wider down- type Higgs profile allows it more overlap with 
the first generation relative to the narrow up-type Higgs profile, which is very small 
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Figure 1: Zero mode profiles for the particular parameter set with l\ = 3.05, I2 = 2.29, 
Z3 = 0.36, lh = —1-7, r u = 1.5, and = 0.67. 

at the distant location of the first family. This results in low energy predictions, 

m M (2 GeV) = 3.6 MeV, m c (m c ) = 1.3 GeV, (2.6) 
m d (2 GeV) = 8.8 MeV, m b (m b ) = 4.3 GeV, 
m s (2 GeV) = 75.0 MeV, m t (m t ) = 166 GeV, 



\V ( 



CKM\ 



( 0.9755 0.221 0.0046 \ 
0.220 0.9748 0.038 
0.0076 0.037 0.9992 



/ 



Comparison of eqs.( |2.4| ) and ( ^.6[ ) indicates that we have succeeded admirably in 
satisfying the low energy flavor measurements. The CKM elements fit comfortably 
into their allowed ranges, and the quark masses are all reasonable. Furthermore, the 
quantities (m u + rrid)/2 = 6 MeV, and m u /md = 0.4 are more rigorously defined 
from Chiral Perturbation Theory (CPT), and fall into the acceptable ranges f22|j . 

This particular model has tan/? = 13. In fact there is some remaining freedom 
to move around in the acceptable ranges for the quark masses and mixings, provided 
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Figure 2: Scatter plot for models with various lh and the resulting tan/?. 

all of the parameters are suitably adjusted together. This can be characterized as 
lh, which shows the greatest freedom to move consistently with the data. In Fig. |2] 
we show a scatter plot of models consistent with low energy data, in the plane of lh 
and the resulting tan f3. As can be seen from the figure, lh can vary between about 
-0.5 and -2.4, with tan/3 running from 45 to about 8 in that range. The cut-off at 
lh = —2.4 occurs because going lower requires r d < 1/2. 

3. Supersymmetry Breaking 

In this section we introduce supersymmetry breaking into the model described 
above. Since our model involves a small extra dimension with SM gauge fields in the 
bulk, it naturally lends itself to the mechanism of gaugino-mediated supersymmetry 
breaking fl"8| , |i~9| . If we use the shining mechanism to break supersymmetry [pOj] , we 
can naturally produce a weak-scale /i term still keeping all coefficients of Oil). We 
shall see that a general consequence of this fi term solution is a non-zero soft mass 
for the down-type Higgs at the high scale. This results in a neutralino LSP in a large 
region of parameter space and thus we avoid bounds on stable charged particles. 

3.1 Gaugino Mediated Supersymmetry Breaking 

It was realized a number of years ago that a very simple set of boundary condi- 
tions on soft terms, namely a non-zero gaugino mass at a high scale, would produce 
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a theory at the weak scale with scalar masses and gaugino masses of the same size, 
suppressed contributions to FCNC and successful EWSB. Models with these bound- 



ary conditions were called "no scale" models [24]. Gaugino mediation [18, 19] gives 
(so far) the simplest realization of this spectrum from a microscopic theory. 

For gMSB to work, we need a small (~ 10 — 100 Planck lengths) extra dimension 
in which the SM gauge fields and their superpartners propagate. We also require 
supersymmetry to break on a 4d hypersurface which is separated in the fifth dimen- 
sion from the hypersurface(s) on which MSSM matter lives. Thus contributions to 
scalar masses via Planck-suppressed operators are also exponentially suppressed as 
Yukawa couplings are in the previous section. 

The effective 4d operator which contributes to gaugino masses is 

d29 Kkr) w " w " (3A) 

where L is the size of the extra dimension, W a is a chiral superfield whose lowest 
component is the gaugino and S is a gauge singlet. The field S lives in the super- 
symmetry breaking sector and has a non-zero VEV in its auxiliary component, F$. 
As a result, this term produces a gaugino mass at the compactification scale of order 

The operators which contribute to squared scalar masses are 



/ 



"mm QtQ * (3 ' 2) 



where Qi are chiral superfields which contain MSSM matter and i,j are flavor in- 
dices. The operators with i ^ j would normally make dangerous contributions to 



FCNC. However, all of the terms in (|3.2|) are multiplied by exponentially suppressed 
coefficients due to the spatial separation between the supersymmetry breaking sec- 
tor and MSSM matter fields, and thus these terms may be ignored provided the 
supersymmetry breaking sector is sufficiently distant from all MSSM matter. This 
suppression was not operative in the gaugino mass operators because the gauginos 
are not localized in y. 

The dominant contribution to scalar masses comes from the renormalization 
group (RG) evolution^ from the compactification scale M c = 1/L to the weak scale. 



The relevant term in the one-loop beta function is 



^Our RG analysis is carried out at two loops with respect to as (with one loop thresholds) and 
at one loop with respect to all other quantities using the beta functions of pfi|. 
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Figure 3: Loop contribution to scalar masses from localized supersymmetry breaking. 



where m 2 ? are the squared scalar masses, are the gaugino masses and Cj(r/) is the 
quadratic casimir for representation r/ of chiral superfield / in gauge group i. As in 
gauge mediation, these contributions are only proportional to gauge couplings, and 
therefor are flavor diagonal and do not contribute to FCNC or CP violation. If we 
take M c ~ Mqut, the loop factor suppression is matched by a large logarithm such 
that scalar masses and gaugino masses are comparable. 

There is also a contribution to the scalar masses which dominates at the com- 
pactification scale. It can be considered a threshold correction from integrating out 
the higher KK modes of the gaugino. It comes from the loop contribution depicted 
in fig. H where gauginos run in the loop and the operators responsible for gaugino 
masses are inserted on the propagators at the location of supersymmetry breaking 
[ Qj , |T9|j . Contributions like this one were calculated in [18[ assuming that the size of 
the extra dimension and the distance between supersymmetry breaking and MSSM 
matter are the same (see also f25|). The latter will not be the case in our model and 
these contributions can be important as we discuss below and in Section |J 

One remaining superpartner mass which we have not yet specified is that of the 
Higgsino. The \x term, which is a superpotential mass term mixing the up- and down- 
type Higgses (H u and H d ) , should be generated dynamically to explain its weak scale 



value required by radiative EWSB. By putting the Higgs fields in the bulk [28, [L£ 



the \i term can be produced via the Giudice-Masiero mechanism |p9|[ . It was also 
pointed out J30| that the shining mechanism [p0|| , which could break supersymmetry 
on a distant hypersurface, could also be used to produce a \i term on our brane. 
Both solutions require somewhat small couplings, the former for soft masses in the 
Higgs sector and the latter for the /i term itself. 



10 



In the standard picture of gMSB, with M c = M GUT and Higgs fields localized 
with MSSM matter, an electrically charged stau is the LSP, a scenario which is 
disfavored by cosmological considerations. The lightest neutralino can become the 
LSP if either the Higgs fields are in the bulk and the soft mass of Hd is somewhat 
larger than the soft mass of H u Jl9| or M c > Mqut by at least an order of magnitude 
so that the scalar masses run significantly above the GUT scale . 

Here we offer modified solutions to the fi and LSP problems which naturally 
fit into our model of flavor. The mechanism allows all couplings to be of 0(1) and 
produces simple high-scale boundary conditions (though the location of supersym- 
metry breaking is somewhat fine-tuned). The weak scale spectrum is in general 
distinguishable from other models of gMSB. 

3.2 Soft Parameters and the fi term 

To break supersymmetry in our model, we will use the shining mechanism of 



Arkani-Hamed, et. al. ||20||. One advantage of shining is that it does not require 



the localization of gauge fields, a difficult task in more than four dimensions [14, 15 



14 


15 


20 


3Cfl 



We discuss two variations of the solution of Schmaltz and Skiba in minimal gaugino 
mediation [[3(J and show in our context how to remove the requirement for a small 



coupling. We also show that in this context that a neutralino can be easily be the 
LSP even if M c = M GUT . 

3.2.1 Scenario 1 

We introduce two chiral superfields, $ and $ c , to the bulk which are singlets under 
the SM gauge groups. These fields together are a hyper mult iplet of the N — 1 super- 
symmetry in 5d, which is conserved by this mechanism up to explicit breaking by a 
source J c which couples to $ in the superpotential. This source is localized at y — l s 
near the MSSM fields. In the absence of further ingredients, the scalar component 
of $ c will acquire a VEV and supersymmetry remains unbroken, supersymmetry is 
broken by introducing another singlet chiral superfield X localized in the bulk far 
from our MSSM matter, and coupling to $ c . Assuming X has a profile i/j^iv) with 
respect to the compact dimension, the Lagrange density can be expressed, 

C = J d A 9 J dy [$(y)&(y) + + |^°(y)f XX^~ 

+ J d 2 6 J dy [^ c (y)(d y + m^y) + J c ^y)S{y-l s )+r, c ^X^ c ] 

+H.c, (3.4) 

where m$, J c , and if are all dimensionful couplings of order the Planck scale to the 
appropriate power. We take tp^iv) normalized such that / dy\ip®(y)\ = 1 in order 
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to produce canonically normalized kinetic terms. The F-term equations are 



If 



X 



\(-d y + m^ c (y) + J c 6(y-l s ) \ 
\(d y + m^(y) + V c ^ x (y)X\ 

dy^Mnv) ■ 



(3.5) 
(3.6) 

(3.7) 



For X localized far from l s , the potential is minimized with a non-zero VEV for <p c : 

(<f) ~ -6{l s - y)J c e- m ^ h - y) . (3.8) 

If tpx(y) is a narrow function localized around the point l x < l s , then the F-term 
conditions (Fj = 0) cannot be simultaneously met and thus supersymmetry is broken. 
The field X now plays the role of the hidden sector singlet" in eqs. ( j3.1| - |3.2D . As 
noted below eq. (|3.2j ), provided X is sufficiently distant from the MSSM matter, it 
will provide negligible soft masses to the sfermions. 

A fi term can be generated [3(1 by adding the following terms to the Lagrangian: 



d 2 6 [J$ c (y)5(y - l x ) + \^° Hu (y)^ Hd (y)H u H d <S>(y)] . 



(3.9) 



Using the notation of Section including the redefinition l { — > and taking 

| l x | ^> 1 we find 

P = j i dyy^X^J^/rJ^esxp [-r 2 u y 2 - r 2 d (y - l h ) 2 - r^y - l x )] 

rlr 2 d l 2 h + r 2 l hr4> + (r^/2) 2 ' 



2r u r d 
r 2 + rl 



cxp 



rl + r 2 . 



M*exp 



r<f>l x - 



r 2 u + r 2 



(3.10) 



where = m^/C,. The first term in the exponential dominates and is negative as 
presumed above. Taking ip x (y) = (2( 2 /tt) 1 ^ 4: e~^ y ~ lx ^ 2 and r x = ( x /(, we get a 
gaugino mass at the compactification scale of 

1/4 



M 1/2 



dy - 



r/ c J c 



7i J ^(M*(M*L) 



; exp [~r 2 x {y - l x f + r^y - l s )] 



(2 



7T 



,1/4. 



T) C J C 



M c exp [r+(l x - l s ) + (r (j> /2r x ) 2 ] . 



exp [-r^lx - l s ) + (r (/> /2r x ) 2 ] 



(3.11) 



"Again the full 5d N = 1 supersymmetry is broken by such couplings. We assume that the theory 
above M* accounts for this specific breaking pattern. A simple high-energy description would be 
to confine X to a brane localized in y, but we shall leave the discussion more general. 
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The last lines of eq.s ( |3.10| ) and (|3.11|) come from taking rf ~ J 2 / 3 ~ ( J c ) 2 / 3 ~ Q x ~ 
M*. By making ~ / s + (a few), the yti term and Mi/ 2 will be the same size. Again, 
a small coupling (required in P0| ) can be generated by an 0(1) distance (in Planck 
units). 

Even though X is fixed far from the MSSM matter fields, loop contributions 
of the type in fig. ^| can still play an important role in the low-energy spectrum 
if the distance between X and the matter fields, Al = (/» — l x ), is a fraction of 
the total size L of the extra dimension. The contribution from these loops were 



calculated in |Tj|] for the case of Al = L. For Al <C L, the loop integrals are cut off 
above the mass of N^k ~ {L/ Al) Kaluza-Klein modes. The integrals require the 5d 
gaugino propagator |18|, [IT] expressed in position space in the fifth dimension and in 
momentum space in the 4 large dimensions, 

where p n = mr/L is the (quantized) momentum flowing the the compact dimension. 
The sum can be carried out and results in a simple expression in terms of hyperbolic 



functions [|yj. Inserting this propagator in the Feynman diagram shown in fig. y the 



contribution to the scalar masses may be estimated as, 

m a.. A ^ r° f catf[q{\ - x)\ 
167r 2 Ji sinh 2 [g] tanh[g] 

where x = Al/L, and we have rescaled q by l/L so that the leading L dependence is 
absorbed into Mi/ 2 and 54, resulting in the 4d effective quantities appearing in the 
equation. 

We have inserted the entire set of SU(5) gauginos into the loops, and assumed 
them all to be massless above the compactification scale. This result should be 
regarded as an order of magnitude estimate, because once the theory begins to look 
5 dimensional, the gauge coupling will experience strong running, and though these 
effects appear in the calculation formally at higher order, it may be important to 
resum them by including the running of the coupling with q. We estimate these 
effects to be a factor of a few. 

These contributions can in principal be larger than one-loop threshold corrections 
and should be included in our high-scale boundary conditions. Distances of Al < L/4 
could lead to contributions which shift the LSP from a stau to a neutralino. What is 
interesting about these contributions is that while they are flavor diagonal they are 
not flavor independent due to the different locations of the generations. For a wide 
range in parameter space (Al vs. L), these contributions are below the bounds on 
additional FCNC and CP violation [II] and leave a distinct imprint on the spectrum. 
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3.2.2 Scenario 2 



For variety, we could move X from where J is to where J c is (near the matter 
fields) and remove J c altogether. This corresponds to exchanging l x and l s above. In 
this way, 0's VEV is responsible for both supersymmetry breaking and the /i term. 
In addition, the coupling 7] C XQ C must be exchanged for rjXQ. In principal this 
would work while giving a non-negligible contribution to Hd from Planck suppressed 
operators. 

For successful gaugino mediation, gravity-mediated contributions to squark and 
slepton masses must be small. Thus, this essentially puts a restriction on how close 
X can be localized to the matter fields. For a squark positioned at Zj, its squared 
mass mj receives a contribution of 



F 2 

m * ~ M* 6XP 



2r 2 



(3.14) 



while an off-diagonal mass squared m?- with squarks localized at U and lj receives a 
contribution 



Ml 



exp 



[{k - Itf + 2rl {{k - l x ) 2 + (/, - l x ) 2 )] 



(3.15) 



which is sufficiently suppressed for distances of order a few units and r x of 0(1). 
Thus, only Hd receives a non-negligible contribution to its soft mass for most of 
parameter space. 

However, the loop contributions described in scenario 1 are too large in this 
scenario because X is so close to the matter fields. This problem can be remedied by 
altering the shining mechanism. If instead of the coupling r]X<& we add the coupling 
\ X X§ 2 then we find for Fx- 



x 



A, 



dyipKyMvf 



(3.16) 



thus requiring l x to have a more intermediate value between 1^ and l s . The \x term 
in this scenario is the same as in equation ( [3.10D if one replaces l x with l s , while the 
universal gaugino mass is 



M 1/2 ~ (2tt) 



1/4 _ 



X X J 2 



y&M*(M*L) 
M c exp [-2r^(l x - l s ) + (r^/r x ) 2 ] . 



exp [-2r^(^ - l 3 ) + (r^/r x ) 2 ] 



(3.17) 



In Section ^ we show that for small enough L, X can give a scalar mass contri- 
bution to Hd which alters the particle spectrum significantly. The loop contributions 
to scalar masses discussed above still play a small but significant role while similar 
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contributions to A terms and B\i are negligible. There are no additional CP violating 
phases coming from soft terms since the phases in M1/2 and \i can be rotated away 
and all of the loop contributions are to diagonal (and therefore real) soft masses** . 



4. Sparticle Spectrum and Phenomenology 

Having introduced the general framework for our model, we now turn to some 
specific numbers. In order to produce a weak scale /i, we fix l s (for a given r s ~ 0(1)), 
such that the scalar VEV of $ reaches the weak scale in the vicinity of the Higgs 
fields. Thus, for a given choice of r s , l s may be determined such that the resulting fi 
term induces the correct EWSB radiatively (though in general one must know what 
the high scale soft masses are in order to know what value of \x that is). 

4.1 Sparticle Spectrum 

Both scenarios for supersymmetry-breaking have similar "gaugino-dominated" 
boundary conditions, in which the soft masses for all scalars are smaller than the 
gaugino masses, with the possible exception of the down-type Higgs, which receives 
the largest gaugino-loop contributions in scenario 1 as well as sizable supergravity 
contributions in scenario 2. 

This hierarchy in the scalar masses at the GUT scale differs in one very im- 
portant way with respect to the standard minimal supergravity inspired (SUGRA) 
models. The fact that we have non-universal scalar masses means that there are 
contributions to the evolution of the soft masses from the U(l) y .D-terms. These 
.D-terms contribute to the beta function of scalar mass mf, 

dt 1 5 167T 2 ' 

S = i m Hu - m H d + Tr [ m Q ~ 2m l + m d ~ m L + m l] ) , I 4 - 1 ) 

where g± is the U(l)y coupling, normalized appropriately for SU(5) unification, Yi 
is (standard) hypercharge for scalar i, and the trace is over the three families. The 
quantity S is an RG invariant, and thus vanishes at all scales if all scalar masses are 
equal at some scale. Thus it contributes nothing to the evolution of SUGRA scalar 
masses. In our model, this term will affect all of the scalar masses, with the effect 
being the most dramatic for the sleptons, which have large hypercharges and receive 
no contributions from the strong coupling. This results in a neutralino LSP in a 
large region of parameter space. 

The condition that Bfi = at the high scale provides us with a particular 
moderate to high value of tan (3 \TJ§ for each choice of M\ji and m^. In fig. ^ 



"This however does not solve the strong CP problem 
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Figure 4: The solid curves are the contours of mH d for a choice of M\j 2 (both at the GUT 
scale) and tan/3. The dashed curves demark regions where t±, cTr, and x? are the LSP. 
The dotted curve marks the value of vn.n d as a function of M^i 2 for the example model of 
scenario 1. 



we show the contours of constant tan/5 in the plane of M 1 / 2 and m Hd , including 
the gaugino-loop contributions given in eq. (|4.2| ). It is remarkable that both the 
mechanism that generates the 4d Yukawa couplings of Section FJ and the scenarios 
for supersymmetry breaking favor the same intermediate to large values of tan (3. 

The resulting LSP in the plane of My 2 an d m Hd is indicated in fig. |] by the 
dashed curves^. From this figure, we see that if M\i 2 is much larger than we 
arrive at a e# LSP because the large M\j 2 results in a relatively large neutralino 
mass, while the selectrons get very little contributions from gaugino loops and the 
D-term contribution to the evolution is small for small m^. On the other hand, 
if m# d is much larger than M\/ 2 , the D-term contribution to the slepton masses 
overpowers the standard gaugino contribution. This results in a T\ LSP (which is 
mostly tl) or a negative m 2 for t% which spontaneously breaks the electromagnetic 



TTThis figure has included scalar masses at the GUT scale introduced in eq. (4.2) below. These 
contributions to the scalar masses at M c do not significantly affect the contours of tan (3, though it 
does generally alter the LSP curves 
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symmetry (for consequences of a generic .D-term, see |j32|). However, for m,H d ~ ^1/2 
the D-term contribution is enough to raise the <3r mass above the lightest neutralino 
mass, and does not push the T\ below it, resulting in a stable xo LSP which provides 
a suitable dark matter candidate, and results in missing transverse energy signatures 
at a hadron collider. 

4.2 Scenario 1 

In this scenario, the gaugino masses arise from a superpotential coupling X $ 
which induces a VEV in Fx- X is located a large distance away from the matter 
fields. Thus, all supergravity contributions to scalar masses are zero because of the 
small overlap between the distant X field and the matter fields. However, there are 
contributions from the gaugino-loops in ( p,13|) which contribute to all of the scalar 



masses, and thus are all proportional to M\j 2 . Thus, all of the boundary conditions 
are effectively specified by My 2 and the position of X. 

As a particular example, for \l x — lh\/L ~ 0.4 one obtains**, 

m\ = 0.014 M? /2 , m 2 Hd = 0.23 M 2 /2 , 
m\ = 0.019 Mh 2 , m\ u = 0.61 M 2 l/2 



f 2 — u ' ul3JH l/2) " l H u — lvl l/2i 

I = 0.049 Ml f2 . 



m \ ~ 0.049 Mf /2 . (4.2) 



Thus, our soft terms for squarks are dominated by the gaugino-mediated contri- 
butions, and are relatively flavor-blind. The off-diagonal squark matrix entries are 
then induced by the CKM rotation from interaction to mass eigenbasis and can be 
estimated for the first two families as, 

^Zp- V us ~ 4 x 10- 4 (4.3) 

This is small enough to avoid the supersymmetry flavor and CP problems generally 
associated with off-diagonal squark masses. The dotted line in fig. |] shows the value 
of mn d for this model. As can be seen, this corresponds to tan/? ~ 10 and has a 
neutralino LSP for M 1/2 < 700 GeV. 

4.3 Scenario 2 

In this scenario gaugino masses are determined by the singlet X which develops 
an auxiliary VEV through superpotential coupling X $ 2 . M\j 2 may be fixed by 
localizing X a suitable distance between the Higgses and the source for $ so that 
Fx/{M 2 L) is a weak scale gaugino mass (we assume that L ~ Z s , or in other words, 



w Though eq.(3.13) is only an estimate of the loop contributions, it preserves the percentage 
difference between scalar masses reasonably accurately. Thus, we present figures to two significant 
digits to better illustrate the mass splittings. 
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that the source for $ is localized on one end of the compact dimension, and the MSSM 
matter and Higgs live roughly at the other end). The X field must be localized far 
enough from the fermions to avoid dangerous supergravity-mediated flavor mixing in 
the soft masses. We can proceed by choosing r s and l s to get a particular weak scale 
fj, and localize X so that an appropriate M1/2 results. For some "typical" numbers 
r s = 2 (and l s = —16.5 so that /i ~ 450 GeV), these two constraints require that 
X be localized in the region l x ~ —4.7 for r x ~ 0.5. This results in a gaugino mass 
of Mi/2 ~ 300, which is the right order of magnitude to provide the correct EWSB 
given the value of \i. 

Computing the supergravity contributions to the soft masses, we find that all of 
the squark and the up-type Higgs receive negligible contributions because they are 
localized too far away from X. On the other hand, the down-type Higgs receives a 
substantial contribution to its soft mass of about rn 2 Hd ~ (250 GeV) 2 because it lies 
relatively close to X. The 5d loop contributions are also generically small because 
they are loop suppressed, and have been chosen for this example to be the same as 
those presented in ( |4.2| ) , and are thus once again safe from the point of view of flavor 
violation. 

By performing a brute force scan through the "reasonable" range of r s and r x 
from 0.5 to 2, we find that this is a generic prediction of the model; we are able to 
obtain any M1/2 and m# d between 100 GeV and 1 TeV, as well as the corresponding 
value of /j, at the high scale required to induce EWSB from these boundary conditions. 
The supergravity contributions to other soft masses are always so small as to be 
negligible, and the gaugino- loop contributions to squark masses are a small factor 
times M\j 2 . 

Having determined that our model easily can realize the plane of M1/2 and 
rriH d from the regions of 100 GeV to 1 TeV, and found that we can accommodate 
the necessary fi for EWSB resulting tan /3, we can switch our discussion from the 
underlying model parameters r s , l s ,r x , and l x and discuss the resulting theory defined 
by a particular choice of M1/2 and at scale M c (which we will take for simplicity 
to be the GUT scale, M GUT ~ 2 x 10 16 GeV). Thus, we may completely specify the 
supersymmetry breaking parameters of scenario 2 by the two free parameters, 

M 1/2 , m Hd (4.4) 

with fi fixed for radiative EWSB and tan /3 determined by the condition that Bfi = 
at M c ~ Mqut- As shown in fig. £|, this results in a particular value of tan/3 
determined by the boundary conditions Bfi = at M c , while \i is determined at the 
weak scale from the observed mass of the Z boson, and depends very strongly on the 
choice of M1/2 and rather weakly on the choice of rriH d - 
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4.4 Collider Signatures 



The resulting weak scale phenomenologies for both supersymmetry breaking sce- 
narios discussed above are similar, and thus we briefly discuss both together here. As 
in all gaugino-dominated models, the sparticle spectrum is mostly dependent on the 
choice of Muz- F° r example, for M1/2 ~ m H d ~ 200 (so tan/3 ~ 12), the resulting 
spectrum has squarks with masses on the order of 430 GeV, charged and neutral 
sleptons with masses between 80 and 150 GeV, gluinos with mass around 500 GeV, 
and weak charginos and neutralinos with masses between 75 and 320 GeV, with the 
lighter states being dominantly gaugino-like. A number of these particles would be 
accessible at Run II of the Tevatron, with a variety of signals (all characterized by 
the missing transverse energy from the Xo LSP). Larger values of Mu% will result in 
a heavier sparticle spectrum, with the qualitative features unchanged; the gluino will 
tend to be the heaviest superparticle, with the squarks somewhat lighter, followed 
by sleptons and weak gauginos. Once the superpartner masses become heavier than 
a few hundred GeV there is insufficient energy to produce them at the Tevatron. 
However, the CERN Large Hadron Collider (LHC) will probe Mu2 up to roughly 
one TeV. Again, this will be realized through a variety of signals, one particular in- 
teresting example of which is the "tri-lepton" signal coming from the leptonic decay 
of the lighter gauginos into the neutralino LSP |33[ . 

As in all models which can be described at low energies by the MSSM, there are 
three neutral Higgs bosons (two CP even and one CP odd) and a pair of charged 
Higgs scalars. The heavy CP even, CP odd, and charged Higgs bosons typically 
have masses that are a few times larger than M1/2 and thus the model exhibits the 
Higgs "decoupling" limit in which the lightest Higgs boson has approximately SM 
couplings. For the moderate values of tan/3 realized by the model, the interesting 
signals of the pseudo-scalar Higgs produced in association with b quarks may be 
observed at Tevatron or LHC p3|| . The lightest Higgs boson typically has a mass in 
the range of 110 to 120 GeV, much of which will be probed by LEP ||35|| , with higher 
masses typically accessible to the Tevatron and/or LHC fl36f. 



5. Conclusions 

We have presented a model in which all masses below the GUT scale are gen- 
erated by localizing fields in a small extra dimension. All dimensionless couplings 
are of order one, all dimensionful couplings are of order the Planck scale and all 
distances are of order the Planck length, save the size of the extra dimension which 
is of order the inverse GUT scale. The quark flavor portion of the model manages 
to beautifully realize the CKM mixings and CP-violating phase observed in nature, 
as well the observed quark masses. It is completely independent of the nature of 
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super symmetry breaking, and thus can be taken as a generic picture of how quark 
masses might arise from a small extra dimension in a supersymmetric context. 

We have explored two different pictures for how the extra dimension might play 
a role in supersymmetry breaking, both of which are variations of gaugino media- 
tion. While successfully incorporating attractive features such as radiative EWSB 
and the possibility of a neutral LSP suitable as a dark matter candidate, they re- 
sult in distinctive boundary conditions for soft masses at the GUT scale, and thus 
result in interesting relations among superparticle masses not seen in other models. 
The separation of the families leaves an imprint on the boundary conditions, pro- 
ducing a striking scalar spectrum distinguishable from other predictive models of 
supersymmetry breaking. 

We have employed a bottom-up approach in constructing the model, taking for 
granted details related to the spontaneous breaking of the GUT symmetry to the SM 
gauge group, specific details concerning the dynamical localization of the fermions, 
and the breakdown of the full N — 1 supersymmetry in the 5d theory to the N = 1 
in the 4d theory. It would be interesting to pursue these last two technical details 
with more rigor. Further, it would be interesting to construct a full GUT theory, to 
explore the possibility that GUT physics stabilizes the size of the extra dimension at 
the GUT scale. Further, the details of the GUT breaking must address the lepton 
masses observed in nature. In fact, our model contains (at least) two gauge singlet 
fields which might play the role of the right-handed neutrino. 

In general, we find that a small extra dimension allows for a vast array of new 
possibilities to solve old problems. While extra dimensions of this size would not 
be accessible at any collider in the near future, their indirect effects on low energy 
phenomena can be dramatic. 
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A. Localizing a Chiral Superfield in an Extra Dimension 

In this appendix we show how a chiral superfield may be localized in an extra 
dimension with an exponential profile. Our starting point is the (on-shell) action for 
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an N = 1 chiral superfield in 5 space-time dimensions, 

+ ^ d 2 £ i tf c [3 B + M(yp + F.c.J | , (A.l) 

where \& and \& c are the left-chiral and charge-conjugated right-chiral iV = 1 4d 
superfield components of the single N = 1 5d chiral superfield d y = d/dy denotes 
the partial derivative with respect to the fifth dimension. This formulation of the 
action is convenient because it is written in the usual N — 1 superfield language and 
thus is manifestly N = 1 supersymmetric. We have written the mass parameter M(y) 
explicitly as a function of the fifth dimension. This could be realized, for example, by 
appropriately coupling an additional chiral superfield to \l/ c ^ and including dynamics 
which give its scalar component a VEV. 



In terms of component fields, eq. (|A.1|) may be written, 

S = J d 4 xdy{F[F L + F* R F R + d^ L d^ L + d^Vfa 

+ iV [7"3„-75 0v-Af(y)]tt 

+ % ( F* R [d y + M{y)\ <\> L + <P R [d y + M(y)\ F L 

+ F* L [d 5 - M(y)\ <p R + 02 [d y - M(y)\ F R )} . (A.2) 

where 4>l, Pl^, F l ((/) r , Pr*&, F r ) are the scalar, spinor, and F component fields of \I> 
(^c), respectively. Expanding these fields in terms of a complete orthonormal basis 
of scalar functions of y, 

P L V(x,y) = J>2(a:)& n (y), P R y(x,y) = J>£(z)/ n (y), 

n n 

<f> L (x,y) = J2^)B n (y), 4> R (x,y) = 

n n 

F L (x, xh) = F ?(x) B n (y)> F R (x, xh) = ^ F n R {x) ^(y), (A.3) 

n n 

we find that the action eq. ( |A.2|) simplifies considerably if one requires, 



f b n (y) ] 1 


' f n (y) 


S»(y) =A n 


T n {y) 


{ B n {y) J 1 


k F n {y) 


f n (y) ) 


\ b n (y) 


F n (y) \ = -\n < 


& n (y) 


T n {y) J 


{ B n {y) 



(A.4) 
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indicating that [b n (y), f n {y)}, [B n (y), T n (y% and [B n (y), F n (y)] are "bosonic" and 
"fermionic" pairs of solutions to a SUSY Quantum Mechanics problem with Q = 
[d y + M(y)] |3^]. These first order differential equations may be combined to give 



second order equations for each individual function, 
[-dl + M\y)-{d y M{y))] 



b n (y) ] 




' b n (y) 


B n {y) 


> = |An| 2 | 


B n {y) 


B n (y) J 






f n (y) y 




\ F(y) 


F n (y) 


| = |A„| 2 < 


F n (y) 






{ f n (y) 



[-d 2 + M 2 (y) + (d y M(y))] I F n (y) } = \\ n \ 2 { F n (y) }. (A.5) 



From this result, we see that the same differential equation determines b n (y), B n (y), 
and B n (y), and another determines f n (y), F n (y), and JF n (?/), and thus there are only 
two independent sets of functions relevant to the four dimensional effective theory 
description (instead of six). This is certainly not surprising; it is an indication that 
4d A^ = 1 supersymmetry remains unbroken. The important point for our purposes 
is that each component field for a given KK mode of the chiral superfield acquires 
the same profile in the extra dimension. Putting this expansion into eq. ( |A.2| ) and 
carrying out the dy integration, we arrive at, 

S = J2 fd 4 x F-* Ft + F%* F£ + d,4>T + 3,<p n R d^l (A.6) 

n J 

+i**kfdp+ \ n ] + i\ n (Fx* $1 + r R * Ft + Ft* r R + <f>? Ft) 

= J2f d4x {J d4 ° (y nj ^ n + y c nj y n c) + d 2 6 i\ n ^ * n + H.c/j | . 

which is a free supersymmetric theory for an infinite tower of massive left- plus right- 
chiral multiplets as well as some number of left- and right-chiral zero mass modes 
(whose number need not be equal J38J). Introducing interactions to this theory 
does not change the end result for the kinetic terms, though generically interactions 
will be induced among all of the various modes, and with varied coupling strengths 
proportional to the overlap of the profiles of all participating modes in the extra 
dimension. 

We will be considering Gaussian profiles for the zero modes, which can be consid- 
ered "generic" in the sense that given any mass function which crosses zero at some 
point in y, for a small region about that point, the function may be approximated as 
linear, M(y) m 2( 2 y, with 2( 2 the slope at the crossing point. An example of such a 
profile is a domain wall arising from a kink soliton. For this choice of M(y), eq. (|A.5| ) 
looks like the Schrodinger equation for a harmonic oscillator with frequency 2( 2 and 
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energy shifted by — ( 2 for the left-chiral and +( 2 for the right-chiral modes. Thus, 
for ( 2 > there is a single left-chiral zero mode with profile, 

b ° {y) = (v) (3)ei " e "' v ' (A ' 7) 

and the higher mode b n (y) are given by the familiar product of exponentials and 
Hermite polynomials. The n > 1 f n (y) functions may be obtained from eq. ( |A.4j) . 
For ( 2 < there is a single right-chiral zero mode. Thus, any crossing of zero in 
M(y) will generally produce a zero mode localized around the zero crossing with 
Gaussian fall-off (the width of the Gaussian being determined by the slope of M(y) 
at the crossing point), provided the deviations from linearity occur sufficiently far 
away from the crossing point. 

eq. (|A.7| ) explicitly includes a phase <p. Such a phase is undetermined by our 
choice of M(y), and does not contribute to the kinetic terms derived in eq. (|A.6|) . 
However, it may play a role in the interaction terms for the 4 dimensional effective 
theory, where careful analysis is required to determine which phases are physical, 
and which may be rotated away by an appropriate redefinition of fields. 
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